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ABSTRACT 

Wc present here the final results of the first spectropolarimetric survey of a small sam- 
ple of active M dwarfs, aimed at providing observational constraints on dynamo action 
on both sides of the full-convection threshold (spectral type M4). Our two previous 
studies (Donati ct al. 2008b; Morin et al. 2008b) were focused on early and mid M 
dwarfs. The present paper examines 11 fully convective late M dwarfs (spectral types 
M5-M8). Tomographic imaging techniques were applied to time-series of circularly 
polarised profiles of 6 stars, in order to infer their large-scale magnetic topologies. For 
3 other stars we could not produce such magnetic maps, because of low variability 
of the Stokes V signatures, but were able to derive some properties of the magnetic 
fields. 

We find 2 distinct categories of magnetic topologies: on the one hand strong ax- 
isymmetric dipolar fields (similar to mid M dwarfs), and on the other hand weak 
fields generally featuring a significant non-axisymmetric component, and sometimes 
a significant toroidal one. Comparison with unsigned magnetic fluxes demonstrates 
that the second category of magnetic fields shows less organization (less energy in the 
large scales), similarly to partly convective early M dwarfs. Stars in both categories 
have similar stellar parameters, our data do not evidence a separation between these 
2 categories in the mass-rotation plane. 

We also report marginal detection of a large-scale magnetic field on the MS star 
VB 10 featuring a significant toroidal axisymmetric component, whereas no field is 
detectable on VB 8 (M7). 
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1 INTRODUCTION 

Magnetic field is a key parameter in theories of stellar for- 
mation and evolution. In cool stars, it powers several ac- 
tivity phenomena, observed on a wide range of wavelengths 
and timescales, which provide a rough proxy of the averaged 
field strength. The well-established rotation-activity relation 
(e.g., Noycs ot al. 1984) supports the idea of a dynamo- 
generated magnetic field in these stars. 
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Larmor (1919) first proposed that the solar magnetic 
field could be induced by plasma motions, already point- 
ing out the importance of shear (fi effect) to generate 
strong toroidal fields and thus explain Hale's polarity law 
of sunspots. I'arkor (1905) completed the basic picture of 
the solar dynamo by introducing the a effect (convection 
made cyclonic by the Coriolis force) to address the issues 
of Cowling's anti-dynamo theorem (Cowling 19.'33) and the 
regeneration of a poloidal field component from a toroidal 
one. The aO. dynamo has thereafter been thoroughly de- 
bated and improved (e.g. Babcock 1961; Lcigliton 1969). 
More recently, helioseismology has been able to probe the 
solar interior and revealed the existence of a thin layer of 
strong shear located at the base of the convection zone: the 
tachocline (e.g., Spiegel & Zahn 1992). Although many as- 
pects of the solar magnetism are still not thoroughly under- 
stood, recent theoretical and numerical studies have pointed 
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out the crucial role of the tachocline as the place of storage 
and amplification of strong toroidal fields (e.g., Ossendrijvcr 
2003; Charbonneau 2005). 

New insight on dynamo processes can be gained from 
the exploration of magnetic fields of cool stars, probing 
the dynamo response in very non-solar regimes of param- 
eters (e.g., fast rotation, deeper or shallower convection 
zone). M dwarfs are particularly interesting since those be- 
low - 0.35 Mq (e.g., CUiabricr t Barafic 1997) are fully 
convective and therefore do not possess a tachocline and 
presumably cannot host a solar- type dynamo. Yet, many M 
dwarfs are known to be active, and these stars follow the 
usual rotation-activity relation (e.g., Dclfossc ct al. 1998; 
Rciucrs & Basri 2007; Kiraga & Stcpicn 2007; West ct al. 
2004). However the strong correlation between X-ray and 
radio luminosities established by Guodcl & Benz (1993) 
for stars of spectral types ranging form F to mid M, is 
no longer valid for very low mass dwarfs which exhibit 
very strong radio emission whereas X-ray emissions dra- 
matically drop (Bergcr 2006). Magnetic fields were also 
directly detected at photospheric level through Zeeman 
effect both in unpolarised (e.g., Saar & Linsky 1985; 
Johns-KruU & Valenti 1996; Rcincrs & Basri 2006) and cir- 
cularly polarised (Donati et al. 2006a) line profiles. 

Spectropolarimetry combined with tomographic imag- 
ing techniques is the optimal technique to investigate the 
magnetic topologies of M dwarfs (see Sec. 3 for more de- 
tails). By recovering the large-scale component of stellar 
magnetic fields, we can provide dynamo theorists with ob- 
servables directly comparable with their modeling (axisym- 
metry, relative importance of the poloidal and toroidal com- 
ponents, characteristic scales...). Previous results have al- 
ready produced strong constraints: Donati ct al. (2006a) 
and Morin ct al. (2008a) (hereafter MOSa) demonstrated 
that the fully convective fast rotator V374 Peg is able to trig- 
ger a strong large-scale axisymmetric poloidal field steady 
on a timescale of 1 year; and exhibits a very low level 

of differential rotation (~ ^^ ). From the analysis of a 
sample of early and mid M dwarfs Donati ct al. (2008b) 
and Morin ct al. (2008b) (hereafter DOS and MOSb respec- 
tively) observed a strong change near the theoretical full- 
convection threshold: while partly convective stars possess 
a weak non-axisymmetric field with a significant toroidal 
component, fully convective ones exhibit strong poloidal ax- 
isymmetric dipole-like topologies. Differential rotation also 
drops by an order of magnitude across the boundary, the 
observed fully convective stars exhibit nearly solid body ro- 
tation. DOS and MOSb also report a sharp transition in the 
rotation-large-scale magnetic field relation close to the full- 
convection boundary, whereas no such gap is visible in the 
rotation-X-ray relation. Considering that X-ray emission is a 
good proxy for the total magnetic energy, it suggests that a 
sharp transition in the characteristic scales of the magnetic 
field occurs near the full-convection limit. This point was 
further confirmed for a few stars by Rcincrs & Basri (2009) 
who report that the ratio of the magnetic fluxes measured 
from circularly-polarised and unpolarised lines draiuatically 
changes across the fully convective limit. Several theoretical 
studies have addressed the challenging issue of dynamo ac- 
tion in fully convective stellar interiors. Durney et al. (1993) 
first proposed that without a tachocline of shear, convec- 



tion and turbulence should play the main role at the ex- 
panse of differential rotation, generating a small-scale field. 
Kiikcr & Riidigcr (1999) and Ghabricr & Kiiker (2006) per- 
formed mean-field modeling of dynamo action in fully con- 
vective stars and found purely non-axisymmetric a^ solu- 
tions, indicating that these objects can sustain large-scale 
magnetic fields. Subsequent direct numerical simulations 
by Doblcr ct al. (2006) and Browning (2008) both realized 
large-scale dynamo action with a significant axisymmetric 
component of the resulting magnetic field. The latter also 
achieved magnetic energy in equipartition with kinetic en- 
ergy and therefore Maxwell stresses strong enough to quench 
differential rotation, resulting in nearly solid-body rotation. 
Despite these recent advances, the precise causes of differ- 
ences between dynamo in fully and partly convective stars 
are not completely understood, and theoretical studies need 
observational guidance. 

In the present paper, we extend our spectropolarimetric 
study to 11 late M dwarfs (spectral types ranging from M5 
to MS). After a brief presentation of the stellar sample and 
of spectropolarimetric observations, we describe the main 
principles of the tomographic imaging process. The Zeeman- 
Doppler Imaging (ZDI hereafter) analysis is then detailed for 
6 stars. For 3 other late M dwarfs, it is not possible to derive 
a definitive magnetic map because of the very low level of 
variability in the Stokes V signatures, but we can still infer 
some information about their magnetic topologies. For the 
very faint stars VB S and VB 10, the noise level is too high to 
allow definite detection of the circularly polarised signatures 
in individual LSD spectra. By averaging the spectra of each 
data set, we marginally detect a large-scale magnetic field on 
VB 10, and show a tentative ZDI reconstruction. We finally 
discuss these results and conclude on the implications of our 
study for the understanding of dynamo processes in fully 
convective stars. 



2 OBSERVATIONS 

2.1 Presentation of the sample 

For this first spectropolarimetric survey, we selected 23 ac- 
tive main-sequence M dwarfs, mostly from the rotation- 
activity study by Dclfossc ct al. (1998), covering a wide 
range of masses and rotation periods (although for a given 
mass the extent in rotation period is rather restricted). In 
the present paper we focus on the low-mass end of the sam- 
ple (O.OS-0.20 Mq): GJ 51, WX UMa, DX Cue, GJ 1245 B, 
GJ 1156 and GJ 3622 are thoroughly studied with tomo- 
graphic imaging techniques. We also present a brief study of 
the large-scale magnetic topologies of GJ 1224, GJ 1154 A, 
GN Leo, VB S and VB 10 (section 10). 

All these stars are known to show signs of activity in 
Ha or X-rays (e.g. Gizis ct al. 2002; Schmitt & Licfkc 2004), 
and photospheric magnetic fields have been measured on 
most of them from the analysis of Zeeman broadening in 
molecular bands (see below). The main properties of the 
sample, inferred from this work or collected from previous 
ones, are shown in Tab. 1. 

Stellar masses are computed from the mass-luminosity 
relation derived by Dclfossc ct al. (2000), based on J band 
absolute magnitude inferred from apparent magnitude mea- 
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Table 1. Fundamental parameters of the stellar sample. Spectral types are taken from Rcid ct al. (1995). Formal error bars derived from 
our study are mentioned between brackets for A/*, Rsini, R^ (Ict) and Prot (3(t), they apply to the last digit of the preceding number. 
See Sec. 2.1 for more details and a discussion about uncertainties. 



Name 


ST 


M^ 


tisini 


Bf 


Prot 


Tc 


Ro 


ijsini 


R* 


i 






(Me) 


(kms-i) 


(kG) 


(d) 


(d) 


(10-2) 


(R©) 


(R©) 


(°) 


GJ 51 


M5 


0.21 (3) 


12 


- 


1.02 (1) 


83 


1.2 


0.24 (2) 


0.22 (3) 


60 


GJ 1156 


M5 


0.14 (1) 


17" 


2.1*' 


0.491 (2) 


94 


0.5 


0.16 (< 1) 


0.16 (1) 


60 


GJ 1245 B 


M5.5 


0.12 (< 1) 


ja 


1.7" 


0.71 (1) 


97 


0.7 


0.10 (2) 


0.14 (< 1) 


40 


WXUMa 


M6 


0.10 (< 1) 


5" 


>3.9'' 


0.78 (2) 


100 


0.8 


0.07 (1) 


0.12 (< 1) 


40 


DXCnc 


M6 


0.10 (1) 


13" 


1.7" 


0.46 (1) 


100 


0.5 


0.07 (1) 


0.11 (< 1) 


60 


GJ 3622 


M6.5 


0.09 (< 1) 


3^ 


- 


1.5(2) 


101 


1.5 


0.09 (6) 


0.11 (< 1) 


60 


GJ 1154 A 


M5 


0.18 (1) 


6'' 


2.1'' 


^ 1.7 


88 


!j 1.9 


- 


0.20 (1) 


- 


GJ 1224 


M4.5 


0.15 (1) 


^3° 


2.7" 


^4.3 


93 


<4.6 


- 


0.17 (1) 


- 


ON Leo 


M5.5 


0.10 (< 1) 


3° 


2.4" 


^ 2.0 


99 


<2.0 


- 


0.12 (< 1) 


- 


VB 8 


M7 


0.09 {< 1) 


5° 


2.3" 


^ 1.0 


101 


;^ 1.0 


- 


0.10 (< 1) 


- 


VB 10 


M8 


0.08 {< 1) 


6" 


1.3" 


^ 0.8 


102 


^0.8 


- 


0.09 (< 1) 


- 



" Rciners & Basri (2007) 

'' Reiners et al. (2009) 

<= Mohanty & Basri (2003) 



surements of 2MASS (Cutri et al. 2003) and Hipparcos par- 
allaxes (ESA 1997). Formal error bars, as derived from 
uncertainties on these measurements are mentioned. The 
intrinsic dispersion of the relation is estimated to be 
lower than 10 %. Radius and bolometric luminosity suited 
to the stellar mass are computed from NextGen models 
(Baraffc ct al. 1998), formal error-bars on stellar mass are 
propagated. The accuracy of these models for active M 
dwarfs is a debated subject (e.g., Ribas 2006), but recent 
studies indicate that the agreement with observations is very 
good for late M dwarfs (Deniory ct al. 2009). For all stars ex- 
cept GJ 51, projected rotational velocities (vsini) are avail- 
able from previous spectroscopic studies. The uncertainties 
on vsini are typically equal to 1 kms~^. For each star, we 
also mention the rotation period (Prot) derived from our 
analysis (see Sec. 3.3), and i?sinj is straightforwardly de- 
duced (with propagated error bar). An estimate of the in- 
clination angle of the rotation axis on the line-of-sight (i) 
is obtained by comparing R sin i and the theoretical radius. 
With this estimate the typical error is of the order of 10° for 
low and moderate inclinations, and 20° for high inclination 
angles, this is precise enough for the imaging process. The 
effect of these uncertainties on the reconstructed magnetic 
maps is discussed in Sec. 3.4. 

We also mention unsigned magnetic fluxes from the lit- 
erature (whenever available in Rciners & Basri 2007 and 
Rciners ct al. 2009) empirically derived from unpolarised 
molecular (FeH) line profiles. These estimates result from 
the comparison with reference spectra of an active and 
an inactive star (corrected for spectral type and rotational 
broadening) which are used to calibrate the relation between 
broadening of magnetically sensitive lines and magnetic flux 
(Reiners & Basri 2006). The authors estimate that the preci- 
sion of the method lies in the 0.5-1 kG range. As this method 
is not sensitive to the vector properties of the magnetic field, 
the Bf values reflect the overall magnetic flux on the sur- 
face of the star. Since Stokes V signatures corresponding to 
neighbouring zones of fields with opposite polarities cancel 
each other, our spectropolarimetric measurements are not 



sensitive to tangled fields and only recover the uncancelled 
magnetic flux corresponding to the large-scale component of 
the magnetic topology. Therefore the ratio of both magnetic 
fluxes is a clear indication of the degree of organization of 
the observed magnetic field. 

The convective Rossby number (Ro) , which is the ratio 
of the rotation period and the convective turnover time, is 
believed to be the relevant parameter to study the impact of 
rotation on dynamo action in cool stars (Noycs ct al. 1984). 
In Table 1, we mention Rossby numbers based on empiri- 
cal convective turnover times derived by Kiraga & Stcpicn 
(2007) from the rotation-activity relation in X-rays. These 
turnover times are ad-hoc fitting parameters that reflect 
more the relation between activity and rotation at a given 
mass than an actual turnover time at a speciflc depth in the 
convection zone. However, the resulting Rossby numbers al- 
low us to compare the effect of rotation on magnetic fleld 
generation in stars having different masses. 



2.2 Instrumental setup and data reduction 

Observations presented here were collected between June 
2006 and July 2009 with the ESPaDOnS spectropolarimeter 
at CFHT . ESPaDOnS provides full coverage of the optical 
domain (370 to 1, 000 nm) in a single exposure, at a resolving 
power of 65, 000, with a peak efficiency of 15% (telescope and 
detector included). 

Data Reduction is carried out with libre-esprit, a 
fully-automated dedicated pipeline provided to ESPaDOnS 
and NARVAL users, that performs optimal extraction of the 
spectra following the procedure described in Donati ct al. 
(1997) that is based on Hornc (1986) and Marsh (1989). 
Each set of 4 individual sub-exposures taken in different 
polarimetric configuration are combined together to pro- 
duce Stokes / (unpolarised intensity) and V (circularly po- 
larised) spectra, so that all spurious polarisation signatures 
are cancelled to first order (Scmcl ct al. 1993; Donati et al. 
1997). In addition, all spectra are automatically corrected 
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for spectral shifts resulting from instrumental effects (e.g., 
mechanical flexures, temperature or pressure variations) us- 
ing telluric lines as a reference. Though not perfect, this 
procedure allows spectra to be secured with a radial veloc- 
ity (RV) internal precision of better than 0.030 kms^^ (e.g., 
Moutou et al. 2007). 

The peak signal-to-noise ratios (S/N) per 2.6 kms^^ 
velocity bin range from 51 to 245, mostly depending on 
the magnitude of the target and the weather conditions. 
An overview of the observations is presented in Table 2, the 
full journal of observations is available in the electronic ver- 
sion. Using the Least-squares deconvolution technique (LSD, 
Donati et al. 1997), polarimetric information is extracted 
from most photospheric atomic lines and gathered into a 
single synthetic profile of central wavelength Ao = 750 nm 
(800 nm for VB 8 and 10). The corresponding effective 
Lande factor geff (computed as a weighted average on avail- 
able lines) is close to 1.2 for all the stars of our sample. The 
line list for LSD was computed from an Atlas9 local ther- 
modynamic equilibrium model (Kurucz 1993) matching the 
properties of our whole sample, and contains about 5, 000 
moderate to strong atomic lines. We notice a multiplex gain 
of about 10 (5 for VB 8 and 10) with respect to the peak 
S/N of the individual spectra of our sample. Although all the 
stars in the sample are active, some exhibit Stokes V LSD 
signatures just above noise level (e.g. DX Cnc), whereas on 
others we detect very strong signatures, with peak-to-peak 
amplitudes as high as 1.8% of the unpolarised continuum 
level (for WX UMa). Temporal variations, due to rotational 
modulation, of the Zeeman signatures is obvious for some 
stars, whereas it is very weak on others, depending e.g., on 
the inclination angle of their rotation axis with respect to 
the line of sight, the complexity and the degree of axisym- 
metry of their magnetic topology. 

For each observation we compute the corresponding lon- 
gitudinal magnetic field (i.e. the line of sight projection) 
from the Stokes I and V LSD profiles through the relation: 



Bi{G) = -2.14 X 10^ 



V V{v) dv 



^0 ffcff c / [Ic- I{v)] dv 



(1) 



(Rees & Semel 1979; Donati et al. 1997; Wade ct al. 2000) 
where v is the radial velocity in the star's rest frame, Ao, 
in nm, is the mean wavelength of the LSD profile, c is the 
velocity of light in vacuum in the same unit as v, Poff is the 
value of the mean Lande factor of the LSD line, and Ic the 
continuum level. 

In the rest of the paper, all data are phased according 
to the following ephemeris: 



HJD = HJDo + ProtE, 



(2) 



where HJDo = 2 453 850 for WX UMa, and HJDo = 
2 453 950 for the other stars; and Prot is the rotational period 
used as an input for ZDI and given in Table 1. 



3 DATA MODELLING 

Zeeman Doppler imaging (Donati & Brown 1997) aims at 
assessing stellar magnetic topologies (at photospheric level), 
from the analysis of time-series of high spectral resolution 



spectropolarimetric observations. In this part we briefly re- 
mind the reader with the main properties of this technique 
and details of our implementation. A more complete descrip- 
tion can be found in M08b and references therein. 

ZDI is an inverse problem, the associated direct problem 
consists in computing the Stokes I and V spectra for a given 
magnetic map. The stellar surface being sampled on a grid of 
~ 1, 000 cells, the local Stokes I and V profiles are computed 
from a model based on Unno-Rachkovsky's equations, for 
a given magnetic map. The magnetic field is decomposed 
into its poloidal and toroidal components and described as a 
set of spherical harmonics-like coefficients, as implemented 
by Donati ct al. (2Q06b). We introduce two filling factors 
that account for subgrid cancellation of polarised signatures 
corresponding to fields of opposite polarities, and allows us 
to accurately reproduce the LSD line profiles (M08b). The 
stellar spectrum is then obtained by a disk integration taking 
into account Doppler shift and limb-darkening. 



3.1 Rotational modulation of polarised lines 

In the presence of a magnetic field one can observe the Zee- 
man effect on spectral lines: (i) unpolarised lines (Stokes I) 
are broadened with respect to the null field configuration 
and (ii) polarised signatures (Stokes Q, U and V) show up. 
Here we only study circular polarisation (Stokes V^), which 
is sensitive to the strength and polarity of the line-of-sight 
projection of the field. Because of the combination of stellar 
rotation and Doppler effect: 

• The contribution of a photospheric region to the stel- 
lar spectrum is correlated with its longitude (at first order). 
Thus, as a magnetic region crosses the stellar disk under 
the effect of rotation, the corresponding polarised signal mi- 
grates from the blue to the red wing of the line. 

• The amplitude of this migration depends on the spot's 
latitude (no migration for a polar spot, maximum migration 
from -V sin i to +v sin i for an equatorial one) . 

• The evolution of the signature during this migration 
(as the angle between the field vector and the line of sight 
changes with rotation) is characteristic of the field orienta- 
tion (e.g. signature of constant polarity for radial field, and 
polarity reversal for azimuthal field). 

Therefore, from a series of polarised spectra providing even 
and dense sampling of stellar rotation, it is possible to re- 
construct a map for the photospheric magnetic field. 



3.2 Magnetic field reconstruction: inverse 
problem 

Starting from a null-field configuration, the series of spec- 
tra computed from a test field is iteratively compared to 
the observed one, until a given Xr level is reached. As the 
problem is partly ill-posed (several magnetic configurations 
can match a data set equally well), the maximum entropy 
solution is selected. The spatial resolution of ZDI depends 
on V sin i as a rule of thumb. The highest degree available in 
the field reconstruction is: 



max( 



27ru sin i 
W 



In 



(3) 



Magnetic topologies of late M dwarfs 5 



Table 2. Synthetic journal of observations. Observation year and number of spectra collected are given in columns 2 and 3. Columns 
4 and 5 respectively list the peak signal to noise ratio (per 2.6 kms~^ velocity bin) and the rms noise level (relative to the unpolarised 
continuum level and per 1.8 kms~^ velocity bin) in the average circular polarisation profile produced by Least-Squares Deconvolution 
(see text) — we precise minimum and maximum values obtained for each observing run. The average value and standard deviation of the 
longitudinal magnetic field (see Eq. 1) and the radial velocity measurements are given in columns 6 and 7. The rotation cycle bounds of 
column 8 are computed with the rotation periods mentioned in Table 1. Complete observation logs are available in the electronic version 
of the article. 



Name 


Year 


'n-obs 


S/N 


f^LSD 
(lO-l/c) 


Be 

(G) 


RV 

(kms-i) 


Cycle 


GJ 51 


2006 


6 


128-165 


7.7-10.1 


-990 (313) 


-5.52 (0.20) 


5.0-9.9 


- 


2007 


9 


159-198 


5.7-7.0 


-1657 (280) 


-6.36 (0.74) 


412.7-418.7 


- 


2008 


9 


118-181 


6.9-10.3 


-1219 (407) 


-6.60 (0.57) 


788.0-794.1 


GJ 1156 


2007 


6 


120-181 


6.8-11.0 


82 (72) 


5.96 (0.24) 


431.9-442.0 


- 


2008 


5 


127-158 


10.4-8.4 


-47 (166) 


5.81 (0.13) 


1091.6-1095.9 


- 


2009 


9 


183-195 


6.5-6.9 


24 (111) 


5.73 (0.24) 


1812.7-1817.0 


GJ 1245 B 


2006 


6 


158-191 


7.1-8.8 


-52 (163) 


5.42(0.11) 


4.2-14.0 


- 


2007 


6 


182-226 


5.7-7.4 


-17 (128) 


5.38 (0.10) 


597.0-601.2 


- 


2008 


10 


138-194 


7.2-10.1 


-6 (67) 


5.46 (0.09) 


1054.6-1057.8 


WXUMa 


2006 


8 


67-142 


19.8-9.6 


-1506 (453) 


70.25 (0.53) 


0.5-4.4 


- 


2007 


6 


115-154 


8.4-11.8 


-1757 (405) 


69.95 (0.06) 


341.8-349.5 


- 


2008 


4 


63-129 


10.4-21.4 


-1811 (271) 


70.15 (0.24) 


755.8-758.6 


- 


2009 


11 


113-163 


8.0-12.5 


-1496 (271) 


69.83 (0.21) 


1214.9-1218.9 


DXCnc 


2007 


5 


119-179 


8.3-12.3 


132 (76) 


10.55 (0.07) 


460.7-471.5 


- 


2008 


7 


90-161 


9.4-16.7 


92 (52) 


10.44 (0.56) 


1160.6-1169.5 


- 


2009 


9 


106-187 


7.8-14.8 


67 (44) 


10.67 (0.10) 


2012.8-2019.6 


GJ 3622 


2008 


8 


128-167 


8.9-11.1 


-32 (29) 


2.27 (0.28) 


397.2-402.6 


- 


2009 


6 


101-162 


9.3-15.0 


-26 (26) 


2.37 (0.05) 


617.3-620.0 


GJ 1154 A 


2007 


6 


118-167 


7.0-10.5 


-714 (76) 


-12.83 (0.11) 


- 




2008 


4 


86-154 


8.3-15.9 


-700 (75) 


-12.92 (0.18) 


- 


GJ 1224 


2008 


12 


51-185 


5.8-19.7 


-563 (41) 


-32.68 (0.04) 


- 


CNLeo 


2008 


4 


172-245 


4.5-6.7 


-691 (54) 


19.62 (0.05) 


- 


VB 8 


2009 


9 


83-107 


15.7-20.0 


29 (53) 


15.39 (0.11) 


- 


VB 10 


2009 


9 


68-80 


22.5-25.7 


58 (61) 


36.23 (0.14) 


- 



where W is the unpolarised local profile width (~ 9 kms"^ 
for an inactive non-rotating M dwarf, M08a). The first term 
in the max function corresponds to the limit of high vsini, 
when line broadening is mostly due to rotation and the 
line profile can actually be seen as one 1-D map of the 
photospheric magnetic field. The second term, ^min is the 
minimum resolution available in the low iisini limit, when 
Doppler shift is small and the information on the field topol- 
ogy mostly comes from the temporal evolution of the shape 
and amplitude of the polarised signatures, ^min can range 
from 4 to 8 mostly depending on the signal to noise ratio of 
the data (DOS, MOSb). 



3.3 Period determination 

Since ZDI is based on the analysis of rotational modulation, 
important inputs of the code are: the projected equatorial 
velocity (vaini), the inclination angle of the rotation axis 
with respect to the line of sight (i) and the rotation period 
(Plot)- As no previous definite period measurement exists 
for any star of our late M subsample, we use ZDI to pro- 
vide a constraint on this parameter. For each star, given the 
II sin i and the theoretical radius (corresponding to the stel- 
lar mass) we can derive an estimate of the maximum value 
for the rotation period as : 



Pmax = 50.6145 X -^, (4) 

V sm I 

where Pmax is expressed in days, 7?* in unit of Rq and 
wsini in kms~^. We test period values in a reasonable range 
(< 1.2 X Pmax), and derive the most probable period as the 
one resulting in the minimum Xr at a given informational 
content (i.e. a given averaged magnetic flux value). We try 
to resolve aliases by comparing the multiple data sets avail- 
able for each star. By fitting a parabola to the resulting Xr 
curve close to the minimum, we derive the optimal value of 
Prot and the associated formal error bar (see MOSb for more 
details). Since several data sets are available for each star, 
we mention the smallest 3cr error bar in Tab. 1. For all the 
studied stars periods inferred from different data sets are 
compatible with each other within the width of the associ- 
ated error bars. 



3.4 Uncertainties on the magnetic maps 

Due to the use of a maximum entropy method, magnetic 
maps reconstructed with ZDI are optimal in the sense that 
any feature present in the map is actually required to fit 
the data. However this method does not allow us to derive 
formal error bars on the reconstructed maps. 

Numerical experiments demonstrate that ZDI provides 
reliable maps and is robust with respect to reasonable un- 
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certainties on various parameters and data incompleteness 
(e.g., Donati & Brown 1997). In addition, our implementa- 
tion based on spherical harmonics and poloidal/toroidal de- 
composition successfully reconstructs global topologies such 
as low-degree multipoles, as well as more complex configura- 
tions (e.g., Donati ct al. 2006b). We also find that ZDI maps 
are robust with respect to the selected entropy weighting 
scheme, provided that phase coverage is complete enough. 

We try to assess the effects of the uncertainties on the 
input parameters (see Sec. 2.1), on the derived magnetic 
quantities in a similar way as Petit et al. (2008). For each 
data set we perform several reconstructions with input pa- 
rameters vsini, i, and Prot varying over the width of the 
error bars and check the resulting map and its properties, 
thus providing a quantitative analysis of the robustness of 
our reconstructions to these uncertainties. We therefore ob- 
tain "variability bars" rather than formal error bars. Given 
the small uncertainties on the rotation period, varying this 
parameter within the 3a" error bar has negligible effect on 
the reconstructed maps. The magnetic quantities listed in 
Tab. 9 to characterize the repartition of magnetic energy 
into different components are affected in different ways by 
variations of the input parameters. In particular, the decom- 
position between poloidal and toroidal energy is very robust, 
the observed variation is less than 10 % of the reconstructed 
magnetic energy. The fraction of magnetic energy in axisym- 
metric modes varies by up to 20 % due to uncertainties on 
input parameters. The most important effect observed is a 
cross-talk between the dipole component and higher degree 
multipoles (in particular the quadrupole), the variation of 
the fraction of energy in the dipole is in the 10-30% range. 
The variability bars on the reconstructed magnetic flux are 
of the order of 30 %. 



4 GJ 51 

We have acquired a total of 24 spectra on GJ 51 split in 
3 series collected on 3 successive years (see Tab. 2). All 
Stokes V spectra exhibit a strong signature of constant po- 
larity (radial field directed toward the star). Temporal vari- 
ation inside each data set is detectable and likely due to 
rotational modulation. It mainly consists of an evolution of 
the signature's amplitude. To our knowledge, no previously 
published measurement of Prot or v sin i exist for this star. 
From our LSD Stokes / profiles we measure mean RV val- 
ues of -5.52, -6.36, and -6.60 kms"^ in 2006, 2007 and 
2008, respectively. These values are in agreement with the 
previously published RV — —7.3 kms~^ (Gizis ct al. 2002). 
Our RV measurements also reveal a drift between the 3 
epochs, too large to be due to convection and that may indi- 
cate a companion orbiting around this star. We also observe 
strong RV temporal variations inside each data set presum- 
ably due to magnetic activity (see Tab. 2), well above the 
intrinsic precision of the instrument (see Sec. 2.2). Our anal- 
ysis of both the unpolarised and polarised spectra leads us 
to wsini = 12 kms~^ and Prot ~ 1.02 d, this is close to 
the 1.06 d period derived from the MEarth^ photometric 
data (J. Irwin, private communication). From these values 



Table 3. Fit achieved for the 3 spectral time-series obtained on 
GJ 51. The imaging process is guided towards a mostly axisym- 
mctric solution, see text. In columns 2 wc give the maximum de- 
gree of spherical harmonics used for ZDI reconstruction. Columns 
3-6 respectively list the initial Xr (i-c without magnetic field), 
the Xr achieved with the imaging process, and the average and 
peak value of the magnetic flux on the reconstructed map. 



Epoch 


^ZDI 


2 
XrO 


Xlf 


{B) 
(kG) 


i^max 

(kG) 


2006 


5 


6.17 


1.00 


1.61 


3.86 


2007 


5 


24.53 


1.00 


1.58 


5.02 


2008 


5 


15.20 


1.00 


1.65 


4.68 



Irwin et al. (2009) 



we infer Psini = 0.24 Rq, whereas for M* = 0.20 Mq evo- 
lutionary models expect P^ = 0.21 Rq. We conclude that 
the inclination angle of the rotation axis with respect to the 
line-of-sight must be high, and set i — 60° for ZDI. We note 
that for Prot ~ 1.02 d none of our data set provide a com- 
plete sampling of the stellar rotation, the best one being the 
2008 data set which covers 30% of the rotation cycle. De- 
spite this poor phase coverage, Stokes V profiles collected at 
the 3 epochs are similar, suggesting that the magnetic field 
is stable and mostly axisymmetric. 

Setting ^zDi ~ 5 for the magnetic field decomposition, 
it is possible to fit each of our 3 data sets with Xr=l- The 
resulting magnetic maps feature a strong non-axisymmetric 
component with a dipole tilted toward the observer, in par- 
ticular those inferred from our 2006 and 2007 spectra. This 
is surprising, it is indeed very unlikely that we observe three 
times GJ 51 at the same phase (when the magnetic pole is 
crossing the line of sight). 

The magnetic map reconstructed by ZDI is highly de- 
pendent on the precise form of the entropy, whereas it is 
generally not the case, indicating that this reconstruction is 
particularly ill-posed. We suggest this is due to the lack of 
information associated with the poor phase coverage, the re- 
constructed solution is strongly infiuenced by the maximum 
entropy constraint: the resulting map is therefore mainly 
composed of a spot of radial field facing the observer. 

We perform another reconstruction of the magnetic field 
of GJ 51, with addition of a priori information in the pro- 
cess, so that it preferentially converges toward a mostly ax- 
isymmetric solution, as far as it allows to fit the data at 
the prescribed x^ level. This is done by putting a strong 
entropy penalty on non-axisymmetric modes, similarly to 
the method used by Donati et al. (2008a) to drive the re- 
constructed topology towards antisymmetry with respect to 
the center of the star. In these conditions, we can also fit 
our 3 data sets with Xr=l- The resulting magnetic maps for 
the 3 epochs are very similar (Fig. 2). The corresponding 
synthetic spectra are plotted along with the data on Fig. 1. 
We find similar results at all epochs: the magnetic topol- 
ogy is almost purely poloidal and axisymmetric, it is mainly 
composed of a very strong dipole aligned with the rotation 
axis (see tables 3 and 9 for more details). The azimuthal and 
meridional component of the field somehow differ between 
the three epochs, but we consider that this is not signifi- 
cant given the weak constraint provided by our data. These 
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Figure 1. Time-scries of Stokes V profiles of GJ 5f , in tiie rest-frame of tfic star, from our 2006 (left-iiand column), 2007 (middle column) 
and 2008 (right-hand column) data sets. Synthetic profiles corresponding to our magnetic models (red lines) are superimposed to the 
observed LSD profiles (black lines). Left to each profile a ztlcr error bar is shown. The rotational phase and cycle of each observation 
is also mentioned right to each profile. Successive profiles are shifted vertically for clarity purposes and the associated reference levels 
{V = 0) are plotted as dotted lines. 



maps are not the only solution permitted by our data, but 
we believe that they represent the most probable one. 

Variability bars are of the same order of magnitude as 
for the other stars. The almost purely poloidal nature of 
the magnetic field is robust to uncertainties on the input 
parameters i and v sin i. In addition, when varying these pa- 
rameters, the topology always features a strong purely ax- 
isymmetric component (i.e. more than 45% of the magnetic 
energy is reconstructed in m = modes) and the main re- 
constructed mode is the radial component of a dipole aligned 
with the rotation axis. However, the values mentioned in 
Tab. 9 for GJ 51 should be considered cautiously as our data 
sets provide a weak constraint and the resulting magnetic 
maps are largely determined by the entropy function used. 
In particular, the high values of longitudinal field measured 
(see Tab. 2) indicate that the large-scale magnetic flux is 
indeed higher than those of mid-M dwarfs studied in M08b. 
But it is not clear whether the large-magnetic flux of GJ 51 
is actually larger than that of WX UMa (see Sec. 7). A defl- 
nite confirmation of the magnetic topology of GJ 51 requires 
multi-site observations to obtain a complete sampling of the 
rotation cycle due to a period close to 1 d. 



5 GJ 1156 

We carried out 3 observing run on the flare star GJ 1156 
— in 2007, 2008 and 2009 — and obtained 20 pairs of 
Stokes / and V spectra. The resulting LSD polarised sig- 
nature is above noise level in nearly all observations. Tem- 



Table 4. Same as Tab. 3 for GJ 1156. 



Epoch tzDI XrO Xrf (-B) B,nax 



(kG) (kG) 



2007 
2008 
2009 



1.95 
2.45 
2.07 



0.95 
1.00 
1.00 



0.06 
0.10 
0.09 



0.32 
0.36 
0.36 



poral variation is obvious on Figure 3, we observe both vari- 
ations of amplitude and polarity. We use vsini — 17 kms~^ 
(Reiners & Basri 2007) which allows us to fit the observed 
polarised and unpolarised profiles, as opposed to the previ- 
ously reported value (usini = 9.2 kms~^, Dclfossc ct al. 
1998). We derive P^ot = 0.491 d, ahhough 1/3 d cannot be 
excluded as a possible period. This is confirmed by MEarth 
photometric periodogram which also exhibits a main peak 
at 0.491 d and another one at 1/3 d (J. Irwin, private com- 
munication). The corresponding Rsini is 0.16 Rq. As evo- 
lutionary models predict R^ = 0.16 R©, we set i=60° for 
ZDI. The rotation period being close to a fraction of day 
our observations (especially those of 2007) do not provide 
an optimal sampling of rotation phases, only our 2008 data 
set provides a reasonable sampling on more than half of the 
rotation cycle. 

The data sets can be fitted down to noise level with 
ZDI for the 3 epochs of observation (see Tab. 4). The cor- 
responding maps of surface magnetic flux are displayed in 
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Figure 2. Surface magnetic flux of GJ 51 as derived from our 2006, (left-iiand column), 2007 (middle column) and 2008 (right-hand 
column) data sets. For GJ 51, the imaging process is adapted to preferentially converge toward a mostly axisymmetric solution in order 
to resolve the ambiguity due to poor phase coverage (see text). The three components of the field in spherical coordinates are displayed 
from top to bottom (flux values labelled in G). The star is shown in flattened polar projection down to latitudes of —30°, with the 
equator depicted as a bold circle and parallels as dashed circles. Radial ticks around each plot indicate phases of observations. 
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Figure 3. Same as Fig. 1 for GJ 1156 2007, 2008, and 2009 data sets (from left to right). 
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Figure 4. Same as Fig. 2 for GJ 1156 2007, 2008, and 2009 data sets (from left to right). 



Figure 4. The 3 maps exhibit similar properties: they mainly 
feature two radial field spots of opposite polarities. The mag- 
netic topologies are thus predominantly poloidal (more than 
80 % of the overall magnetic energy in poloidal modes at all 
epochs) but feature a significant toroidal component (in par- 
ticular in 2007 and 2008), and non-axisymmetric (more than 
80 % of the overall magnetic energy in modes with azimuthal 
number m > £/2) as was expected from the polarity reversal 
of the V signature inside each data set. The lower average 
magnetic flux, as well as the weakness of the spot of negative 
polarity (incoming field lines, in blue) , on the 2007 map can 
be attributed to poor phase coverage. 

The reconstructed maps are quite robust to uncertain- 
ties on the input parameters (see Sec. 3.4). In particular, 
when varying these parameters, the fraction of magnetic en- 
ergy reconstructed in non-axisymmetric modes {m > 1/2) 
is always higher than 70%. Using the alternative values 
frot=0.33 d and i=40°, the reconstructed magnetic maps 
do not change significantly. All the quantities listed in ta- 
ble 9 vary by less than 10 % of the total magnetic energy. 
The reconstructed magnetic fiux variations range from 10 to 
15 %. Considering these values therefore does not affect our 
conclusions. 



6 GJ 1245 B 

The M5.5 dwarf GJ 1245 B was observed during 3 succes- 
sive years, for a total of 22 pairs of Stokes / and V spectra. 
The circularly polarised LSD signatures (see Fig. 5) have 
a moderate amplitude and exhibit strong variability (am- 
plitude, shape and polarity) during an observation run — 
presumably due to rotational modulation. Variability also 



Table 5. Same as Tab. 3 for GJ 1245 B. 



Epoch 


izDi 


9 
XrO 


xh 


(kG) 


J^max 

(kG) 


2006 


4 


4.41 


1.00 


0.17 


0.47 


2007 


4 


4.92 


1.10 


0.18 


0.58 


2008 


4 


1.81 


1.00 


0.06 


0.22 



seems important between the different epochs, in particular 
in the 2009 data set the average amplitude of the signatures 
is significantly lower than at previous epochs, this is also 
visible in the longitudinal field measurements (see Tab. 2). 

We measure a mean RV = 5.4 kms~^ with typical 
dispersions of 0.1 kms~^ (see Tab. 2). This is compatible 
with the mean value of 5 kms~^ reported by Dclfossc ct al. 
(1998). We use vsini = 7 kms"^ (Dclfossc ct al. 1998 ; 
Rcincrs & Basri 2007); and find a rotation period of 0.71 d 
corresponding to a peak in the photometric periodogram 
produced by the HATNet ~ survey (J. Hartman, private 
communication). With these values of wsini and period, we 
find it sin i = 0.10 R©, as the NextGen evolutionary model 
predicts Ri, = 0.14 Rq we set i = 40° for our study. 

Running ZDI on the LSD time-series, with the afore- 
mentioned parameters, we can achieve a good fit for the 3 
epochs (see Fig. 5 and Tab. 5). The reconstructed magnetic 
field significantly evolves between two successive epochs, in 
particular the reconstructed magnetic flux has strongly de- 
creased between our 2008 and 2009 observations. However 

2 Bakos et al. (2004) 
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Figure 5. Same as Fig. 1 for GJ 1245 B 2006, 2007, and 2008 data sets (from left to rigiit). 
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Figure 6. Same as Fig. 2 for GJ 1245 B 2006, 2007, and 2008 data sets (from left to rigiit). 
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Table 6. Same as Tab. 3 for WX UMa. 



Epoch IzDl 



XrO 



X^rf 



\-D/ J^max 

(kG) (kG) 



2006 


4 


8.98 


0.90 


0.89 


3.82 


2007 


4 


20.48 


1.05 


0.94 


4.88 


2008 


4 


12.42 


1.00 


1.03 


4.55 


2009 


4 


14.58 


1.15 


1.06 


4.53 



Table 7. Same as Tab. 3 for DX Cnc. 



Epoch 


izDi 


2 
XrO 


xh 


{B) 
(kG) 


J^max 

(kG) 


2007 


6 


1.70 


1.00 


0.11 


0.22 


2008 


6 


1.42 


1.00 


0.08 


0.20 


2009 


6 


1.37 


1.00 


0.08 


0.18 



the magnetic topologies feature similar properties: strong 
spots of radial field (although more than 40 % of the mag- 
netic energy lies in non-radial field structures); a mostly non- 
axisymmetric field (more than 50% of the magnetic energy) 
and a significant toroidal component (between 15 and 20% 
of the total magnetic energy). The mainly non-axisymmetric 
nature of the magnetic field in 2006 and 2008, as well as the 
presence of a significant toroidal component at all epochs are 
robust to uncertainties on stellar parameters (see Sec. 3.4). 



field are very weak. The evolution of the magnetic field over 
successive years is very weak, the maps are strikingly simi- 
lar. These conclusions are robust to uncertainties on stellar 
parameters (in particular i and wsini). When varying these 
parameters over the width of their respective error bars (see 
section 3.4) the reconstructed magnetic field of WX UMa 
is always almost purely poloidal, mostly axisymmetric, and 
the main mode is the radial component of a dipole aligned 
with the rotation axis. 



7 WX UMA=GJ 412 B 

Between 2006 and 2009, we observed the M6 dwarf WX UMa 
during 4 runs, collecting a total of 29 spectra. The LSD 
Stokes V profiles are very similar throughout the data set: 
a very strong (the peak-to-peak amplitude is close to 2% 
of the unpolarised continuum level) simple two-lobbed sig- 
nature of negative polarity (i.e. corresponding to a lon- 
gitudinal field directed toward the star). Temporal evo- 
lution inside each data set, presumably due to rotational 
modulation, is noticeable though weaker than what we ob- 
serve on GJ 51. We measure average values ranging from 
RV = 69.95 to 70.25 kms^^ with a jitter ranging from 
0.06 to 0.53 km s~^ (depending on the observation epoch, 
see Tab. 2). This is in agreement with RV = 68.886 kms~^ 
(single precise measurement by Nidcvcr ct al. 2002, on 
GJ 412A). We use vsmi = 5.0 kms'^ (Rcincrs ct al. 2009) 
which accounts for Zeeman broadening and results in bet- 
ter agreement with our Stokes / and V LSD profiles than 
the previous value of vsini = 7.7 kms~^ inferred from cor- 
relation profiles (Delfosse et al. 1998). From the circularly 
polarised LSD profiles we infer Prot = 0.74 d. With this ro- 
tation period, our 2007 data cover half of the rotation cycle, 
and the 3 other data sets result in a good phase coverage. 
The comparison of the resulting i?sini — 0.073 R© with 
R* = 0.12 R© (from theoretical models) indicates an inter- 
mediate inclination angle, we set i = 40°. 

Although only the 2006 data set can be fitted below 
Xr = 1-0 (see Tab. 6), the ZDI synthetic Stokes V pro- 
files match well the evolution of the LSD signatures for all 
epochs, as shown in Fig. 7. The corresponding magnetic 
maps are presented in Fig. 8: they all feature a strong polar 
cap of radial field (of negative polarity, i.e. field lines directed 
toward the star) reaching a maximum fiux of approximately 
4 kG, whereas the magnetic fiux averaged over the visible 
fraction of the star is about 1 kG. Azimuthal and meridional 
field structure are much weaker. The topology is very simple, 
modes with degree f > 4 can be neglected, the dipole modes 
encompass more than 60% of the reconstructed energy at all 
epochs. Toroidal and non-axisymmetric components of the 



8 DX CNC=GJ 1111 

We carried out 3 observation runs on DX Cnc between 2007 
and 2009, resulting in 21 pairs of Stokes / and V spec- 
tra. The LSD polarised profiles are displayed in Fig. 9, 
the Zeeman signatures have low amplitudes but are defi- 
nitely detected in several spectra. The amplitude and shape 
of the circularly polarised line dramatically evolve during 
each observing run, presumably due to rotational modula- 
tion. From the LSD profiles, we measure average RV rang- 
ing from 10.44 to 10.67 kms^^, with a jitter strongly de- 
pending on the observing epoch (see Tab. 2), in agreement 
with the previously published values (RV = 9 kms^^ in 
Delfosse ct al. 1998; RV = 10.1 kms"^ in Mohanty & Basri 
2003). We use vsini = 13 kms"^ (Rciners & Basri 2007), 
and Prot = 0.46 d (inferred from our data). As the resulting 
Psini = 0.12 R© is already higher than P* = 0.11 R© pre- 
dicted by theoretical models, we assume a high inclination 
angle of the rotation axis and set i = 60° for the imaging 
process. Our 2007 and 2008 data sets provide a reasonable 
phase coverage and the 2009 one results in a good sampling 
of the stellar rotation. 

Both the mean longitudinal field and the standard de- 
viation from this value get weaker from one observing run 
to the next one, indicating intrinsic variability of the mag- 
netic field. Using the ZDI tomographic imaging code, we can 
fit the 3 data sets down to Xr = 1-0, the resulting magnetic 
fields are presented in Fig. 10. Although for the three epochs, 
the reconstructed magnetic topologies feature a significant 
non-axisymmetric component, they significantly differ from 
each other: (i) the fraction of magnetic energy reconstructed 
in the toroidal component grows from 7% in 2007 to 38% in 
2009; (ii) the two main spots of radial magnetic field seem 
to evolve between 2007 and 2008, and finally in 2009 only 
one region of strong radial field remains, (iii) The averaged 
magnetic fiux decreases from 110 G in 2007 to 80 G in 2008 
and 2009. This last point is strengthened by the fact that 
for a given topology, ZDI recovers less magnetic flux for a 
dataset providing partial phase coverage (due to the max- 
imum entropy constraint), whereas here the larger fiux is 
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Figure 8. Same as Fig. 2 for WX UMa 2006, 2007, 2008 and 2009 data sets (from left to right). 
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Figure 9. Same as Fig. 1 for DX Cnc 2007, 2008, and 2009 data sets (from left to right). 






■u-ViL-inai ":td 




Wtrifl fl'^ : >I4 










Figure 10. Same as Fig. 2 for DX Cnc 2007, 2008, and 2009 data sets (from left to right). 
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Table 8. Same as Tab. 3 for GJ 3622. 



Epoch 






(kG) (kG) 



2008 
2009 



2.08 0.95 0.05 0.11 
1.68 0.95 0.06 0.11 



inferred from the dataset providing the poorest sampling of 
stellar rotation. The presence of a strong toroidal compo- 
nent on DX Cnc in 2008 and 2009 is a robust result. In 
particular, when varying the input parameters within their 
respective error bars (see Sec. 3.4) the toroidal component 
always accounts for at least 30 % of the reconstructed mag- 
netic energy. 



9 GJ 3622 

GJ 3622 was observed in 2008 and 2009, we collected 8 
and 5 sequences, respectively. In spite of relatively low S/N 
(due to a low intrinsic luminosity), Stokes V signatures are 
clearly detected in several spectra and variations are un- 
doubtedly noticeable (see Fig. 11). From our data sets we 
infer Prot = 1-5 d, and use vsini = 3 kms^^ reported by 
Mohanty & Basri (2003), resulting in Rsini = 0.09 Rq. Our 
data sets provide a reasonable phase coverage with this pe- 
riod. As evolutionary models predict i?* = 0.11 Rq for a 
0.09 Mq M dwarf, we set the inclination angle to 60° for the 
imaging process. 

Setting (.max ~ 2, given the very weak signal detected, 
it is possible to fit our 2 data sets down to noise level (see 
Tab. 8). The corresponding magnetic maps (see Fig. 11) are 
very similar. A radial field spot of negative polarity, i.e. field 
lines directed toward the star, is located at mid-latitudes. 
Magnetic flux reaches up to 110 G in this region. Azimuthal 
and meridional fields are much weaker, the toroidal com- 
ponent represents less than 10 % of the overall magnetic 
energy at both epochs. The reconstructed topology is close 
to a tilted dipole: less than 10 % of the magnetic energy is 
reconstructed in ^ = 2 modes, and the axisymmetric compo- 
nent stands for more than 80 % of the energy content. The 
data being very noisy the simple topology reconstructed by 
ZDI likely reflects the lack of information in the polarised 
spectra. 



10 OTHER STARS 

For five stars of our sample we collected time series of po- 
larised spectra but could not produce a definitive magnetic 
map. For GJ 1154 A, GJ 1224 and CN Leo (GJ 406) we de- 
tect very strong and simple signatures (see Fig. 12). To our 
knowledge, no rotation periods have been measured for these 
stars and our data sets do not allow us to conclude, either be- 
cause of low intrinsic variability of the Zeeman signature or 
poor phase sampling. Although we cannot compute a mag- 
netic map for these stars, the collected spectra unmistakably 
show that they host very strong large-scale magnetic fields 
(longitudinal fields are about 600 G). Strong magnetic fields 



(total magnetic fluxes in the 2-3 kG range) have been pre- 
viously detected on these stars by Rcincrs & Basri (2007) 
and Rxnncrs et al. (2009) from the analysis of unpolarised 
spectra (see Tab. 1). The simple signatures (two-lobbed anti- 
symmetric) featuring very low variability also clearly suggest 
that these fields are mostly poloidal, strongly axisymmetric 
and presumably dominated by low degree modes, similar to 
what we observe on GJ 51 and WX UMa for instance. The 
low dispersions of longitudinal fields and RV values in each 
data set may indicate that phase sampling is loose and thus 
rotation period close to a fraction of day, or/and that these 
stars are observed nearly pole-on. 

For the faintest stars of the sample VB 8 (GJ 644 C) 
and VB 10 (GJ 752 B), the Stokes V signatures are too 
weak to be definitely detected in individual LSD spectra. 
The initial Xro ^^^ respectively equal to 1.089 and 1.150. 
LSD signatures observed on these stars are shown on Fig- 
ure 13. By averaging all the LSD profiles of a data set, the 
noise level is decreased but only features visible on all spec- 
tra — corresponding to the axisymmetric component of the 
field, if rotation sampling is even — remain visible. The re- 
sulting signal shown in Fig. 13 (bold red line) is processed 
with a zero phase shift low-pass filter to remove the fre- 
quencies higher than permitted by the instrumental profile 
(width of 4.6 kms"\ or 2.5 LSD pixels). For VB 8, the av- 
eraged LSD profile does not feature any significant signal 
(Xr=0.99), indicating that the axisymmetric component of 
the magnetic field of this star is too weak to be detected. 
The averaged profile of VB 10 features a weak but distin- 
guishable signature corresponding to Xr ~ 1-90. It suggests 
the presence of a large-scale magnetic field having a signifi- 
cant axisymmetric component, but further observations are 
needed to confirm this point. From observations in unpo- 
larised light, Reinors A; Basri (2007) report total magnetic 
fluxes of 2.3 and 1.3 kG on VB 8 and VB 10, respectively 
(see Tab. 1). The very weak Stokes V signatures observed 
here (with corresponding maximum longitudinal flelds of the 
order of 100 G) suggest that the magnetic fleld of these stars 
is mainly structured on small spatial scales. 

Performing a ZDI analysis on the VB 10 time series, we 
find 2 possible rotation periods: 0.52, and 0.69 d, the sec- 
ond being favored by photometric measurements (MEarth 
project, J. Irwin, private communication). Figure 14 shows 
the fit achieved for Prot = 0.69 d, i = 60 ° and the 
corresponding magnetic maps. As expected from the sig- 
nature shape (non antisymmetric with respect to the line 
centre), the reconstructed magnetic field exhibits a signifi- 
cant axisymmetric toroidal component. A non-axisymmetric 
poloidal field (tilted quadrupole, mode 022) is also recon- 
structed to fit the Stokes V component that varies along 
rotation. 



11 DISCUSSION AND CONCLUSION 

We present the final part of our exploratory spectropolari- 
metric survey of M dwarfs, following D08 and M08b (respec- 
tively concentrating on mid and early M dwarfs) we focus 
here on the low mass end of our sample. For 6 stars, it is 
possible to apply ZDI techniques to the time-series of cir- 
cularly polarised spectra and thus to infer the large-scale 
component of their magnetic topologies. The properties of 
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Figure 11. 2008 and 2009 spectra and magnetic maps of GJ 3622. See Fig. 1 and 2 for more details. 
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Figure 12. Stolccs / (lower panels) and V (upper panels) LSD signatures of GJ 1154 A, GJ 1224 and ON Leo, from left to right. In each 
panel all the profiles of the time series arc plotted as superimposed grey lines, and the average profile is shown in red. Vertical dotted 
lines represent the line center (in bold) and the approximate limits of the line. In Stokes V plots ztlcr levels (corresponding the individual 
spectra) are shown as dashed lines, and the reference level as a dotted line. 



the reconstructed topologies of these stars are presented in 
Table 9. For the remaining 5 objects, the data sets do not 
permit such a study, it is however possible to retrieve some 
constraints about their magnetic properties. 

Two stars of the subsample (namely GJ 51 and 
WX UMa) exhibit large-scale magnetic fields very similar to 
those observed by M08b on mid M dwarfs, i.e. very strong, 
axisymmetric poloidal and nearly dipolar fields, with very 
little temporal variations. For three stars for which we can- 
not perform a definitive ZDI reconstruction (GJ 1154 A, 
GJ 1224 and CN Leo) our data strongly suggest similar 
topologies. From the observations of WX UMa, we demon- 



strate that the timescale of temporal evolution in the mag- 
netic topologies of these stars can be larger than 3 years, 
whereas previous observations by M08a and M08b were 
based on observations spanning only 1 year. 

The other stars for which we can reconstruct the large- 
scale magnetic topologies, are clearly different. They are 
weaker than those of the first category, and generally feature 
a significant non-axisymmetric component, plus a significant 
toroidal component (although the field is always predomi- 
nantly poloidal). Temporal variability is also noticeable, in 
particular for GJ 1245 B our data indicate unambiguously 
that the magnetic field strength has dramatically decreased 
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between our 2007 and 2008 observations. Our conclusions 
are robust to uncertainties on stellar parameters (Prot , i and 
Dsini). When varying these parameters over the width of 
their respective error bars (see section 3.4), the main prop- 
erties of the reconstructed magnetic topologies remain un- 
changed. 

These results are presented in a more visual way in 
Fig. 15. Previous studies by DOS and M08b have revealed 



strong evidence that a clear transition occurs at approxi- 
mately 0.5 M0 i.e. more or less coincident with the transi- 
tion to a fully convective internal structure. The situation 
here is different, we find stars with similar stellar parame- 
ters that exhibit radically different magnetic topologies. On 
Fig. 15, WX UMa is the only star below 0.2 M© to host 
a mid-M-dwarf-Uke field. Whereas DX Cnc and GJ 1245 B 
are very close to it in the mass-rotation plane they feature 
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Table 9. Magnetic quantities derived from our study. For each star, different observation epochs are presented separately. In columns 
2-4 we report quantities from Table 1, respectively the stellar mass, the rotation period, and the effective Rossby number. Columns 5, 
6 and 7 mention the Stokes V filling factor, the reconstructed magnetic energy and the average magnetic flux. Columns 8-fl list the 
percentage of reconstructed magnetic energy respectively lying in poloidal, dipole (poloidal and I = 1), quadrupole (poloidal and £ = 2), 
and octupolc (poloidal and £ = 3) modes. In column f2, we mention the percentage of magnetic energy reconstructed in axisymmctric 
modes (defined as m < £/2) and the percentage of poloidal energy in axisymmctric modes. See section 3.4 for a discussion on the 
robustness of magnetic map reconstruction and the uncertainties associated with the derived quantities. 



Name 


Mass 


Prot 


Ro 


fv 


(B^) 


{B) 


pol. 


dipole 


quad. 


oct. 


axisymm. 




(M©) 


(d) 


(10-2) 




(105 G2) 


(kC) 


{%) 


(%) 


(%) 


(%) 


(%) 


GJ 5fi (06) 


0.21 


1.02 


1.2 


0.12 


38.6 


1.61 


99 


96 





2 


91/91 


(07) 


- 


- 


- 


0.12 


31.3 


1.58 


99 


92 





6 


77/77 


(08) 


- 


- 


- 


0.12 


32.6 


1.65 


97 


92 


1 


3 


89/89 


GJ lf562 (07) 


0.14 


0.49 


0.5 


1.0 


0.06 


0.05 


88 


30 


26 


19 


6/3 


(08) 


- 


- 


- 


1.0 


0.19 


0.11 


83 


41 


28 


11 


20/12 


(09) 


- 


- 


- 


1.0 


0.13 


0.10 


94 


54 


24 


10 


2/1 


GJ 1245 B (06) 


0.12 


0.71 


0.7 


0.06 


0.44 


0.17 


80 


45 


14 


13 


15/9 


(07) 


- 


- 


- 


0.10 


0.49 


0.18 


84 


46 


27 


7 


52/53 


(08) 


- 


- 


- 


0.10 


0.06 


0.06 


85 


33 


25 


19 


20/18 


WX UMa (06) 


0.10 


0.78 


0.8 


0.12 


16.08 


0.89 


97 


66 


21 


6 


92/92 


(07) 


- 


- 


- 


0.12 


24.42 


0.94 


97 


71 


13 


3 


92/94 


(08) 


- 


- 


- 


0.12 


23.53 


1.03 


97 


69 


19 


6 


83/85 


(09) 


- 


- 


- 


0.12 


37.54 


1.06 


96 


89 


2 


2 


95/96 


DX Cnc (07) 


0.10 


0.46 


0.5 


0.20 


0.17 


0.11 


93 


69 


11 


9 


77/77 


(08) 


- 


- 


- 


0.20 


0.09 


0.08 


73 


31 


25 


10 


49/34 


(09) 


- 


- 


- 


0.20 


0.09 


0.08 


62 


42 


11 


4 


70/61 


GJ 3622 (08) 


0.09 


1.5 


1.5 


1.0 


0.04 


0.05 


96 


90 


7 


- 


73/72 


(09) 


- 


- 


- 


1.0 


0.05 


0.06 


93 


84 


9 


- 


80/78 



^ For GJ 51, the imaging process is weakly constrained by our data sets due to poor phase coverage (see Sec. 4). 

^ For GJ 1156 the alternative rotation period 0.33 d cannot be definitely excluded (see Sec. 5). In this case, the reconstructed 

topologies remain similar for the 3 epochs, and our conclusions arc not affected. 



fields with very diflferent properties. This observation may 
be explained in several ways. For instance, another parame- 
ter than mass and rotation period, such as stellar age, may 
play a role. In our sample we indeed notice that most stars 
below 0.15 M© that exhibit a weak complex field belong to 
a young kinematic population according to Dclfossc et al. 
1998 (GJ 1156, DX Cnc and GJ 3622), whereas those hosting 
a strong dipolar field (WX UMa, GJ 1224 and CN Leo) be- 
long to older kinematic populations (old disk and old/young 
disk). This hypothesis requires further investigation. One 
could also imagine, for instance, that the magnetic fields of 
very low mass stars may switch between two different states 
over time. This hypothesis is supported by the fact that 
for one of the stars in the weak and complex field regime 
(GJ 1245 B) we observe a dramatic variation of the mag- 
netic flux on a timescale of one year which may indicate that 
the magnetic field of these objects may go through chaotic 
variations and eventually switch between the two categories 
of field actually observed. However no such switch has been 
observed in our sample. We observe 5 objects in the strong 
dipole field category (GJ 51, GJ 1154 A, GJ 1224, CN Leo 
and WX UMa) and 6 in the weak field category (GJ 1156, 
GJ 1245 B, DX Cnc, GJ 3622, VB 8 and VB 10), indicat- 
ing that stars would spend as much time in both states. No 
star is observed in an intermediate state, suggesting that a 
putative transition would be fast. This hypothesis may be 
investigated through the analysis of long-term radio moni- 
toring, since radio emission would be presumably strongly 



impacted by such a dramatic change of the stellar magnetic 
field. Recent observations of ultracool dwarfs reveal a long- 
term variability of the activity indices (e.g., Antonova ct al. 
2007; Bcrgcr ct al. 2010) that may support this view. 

On the left panel of Fig. 16, we plot the reconstructed 
magnetic flux as a function of the Rossby number. As men- 
tioned in DOS, stars with M* > 0.4 M© follow the ex- 
pected rotation-magnetic field connection, with saturation 
for Ro < 0.2. Whereas for Af* > 0.4 M©, we only observe 
objects in the saturated regime, with a significantly stronger 
saturation magnetic flux. All the stars of the late M subsam- 
ple have Rossby numbers below 2 x 10^^ and are thus ex- 
pected to be in the saturated regime. As in Fig. 15 the stars 
studied here can be divided into two distinct categories. The 
flrst one is composed of stars hosting a very strong magnetic 
fleld that lie in the saturated part of the rotation-activity re- 
lation, similar to mid M dwarfs. We notice that significantly 
higher magnetic fluxes are reconstructed for WX UMa and 
GJ 51 than for mid M dwarfs studied by M08b, although 
all these stars seem to lie in the saturated dynamo regime. 
For the second category of late M dwarfs (having a weak 
complex fleld), we recover less magnetic flux (even less than 
for some early M dwarfs studied by DOS), in spite of their 
very low Rossby number. Super-saturation is unlikely to be 
the explanation since stars of both categories have similar 
Rossby numbers. As mentioned for Fig. 15, we observe no 
object in an intermediate state. The aforementioned vari- 
ability of GJ 1245 B is also prominent in this plot. 
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Figure 15. Properties of the magnetic topologies of our sample of M dwarfs as a function of rotation period and stellar mass. Larger 
symbols indicate larger magnetic fields while symbol shapes depict the different degrees of axisymmetry of the reconstructed magnetic field 
(from decagons for purely axisy mmetric fields to sharp stars for purely non axisymmetric fields) . Colours illustrate the field configuration 
(dark blue for purely toroidal fields, dark red for purely poloidal fields and intermediate colours for intermediate configurations). Solid 
lines represent contours of constant Rossby number Ro = 0.1 and 0.01 respectively corresponding approximately to the saturation and 
super-saturation thresholds (e.g., Pizzolato ct al. 200.3). The theoretical full-convection limit (M* ~ O.SSMq, Cliabricr & Baraffc 1997) 
is plotted as a horizontal dashed line, and the approximate limits of the three stellar groups discussed in the text are represented as 
horizontal solid lines. Stars with M* > 0.45 M0 are from DOS, whereas those with 0.25 < M* < 0.45 M0 are from M08b. For GJ 1245 B 
symbols corresponding to 2007 and 2008 data sets are superimposed in order to emphasize the variability of this object. Uncertainties 
associated with the plotted magnetic quantities are discussed in section 3.4. 



On the right panel of Fig. 16, magnetic flux inferred 
from Stokes V measurements are compared to those derived 
from Stokes /. Again, stars with masses below and above 
0.4 Mq studied by DOS and M08b clearly form two separate 
groups. In partly-convective stars only a few percents of the 
magnetic flux measured in / is detectable in V similarly to 
what is observed for the Sun, while this ratio is close to 15% 
in fully convective ones. This indicates a higher degree of 
organization of the fleld in fully convective mid M dwarfs, 
with more magnetic flux in the spherical harmonics of lowest 
degree. The ratios of about 30% plotted for WX UMa are 
upper limits (since the flux based on Stokes 7 is a lower 
limit), it is thus not clear if this star differs from mid M 
dwarfs in this respect. This high value however indicates a 



high degree of organization. For DX Cnc, GJ 1245 B and 
GJ 1156, the ratio of Stokes V and / fluxes is closer to the 
early M dwarfs value. Therefore the magnetic fields of the 
weak field category of late M dwarfs share similar properties 
with that of partly convective M dwarfs. 

These results on the magnetic topologies of M dwarfs 
also suggest that dynamo processes in low-mass main se- 
quence stars and pre-main sequence stars may be similar. 
Indeed the first spectropolarimetric results on young stars 
show that the fully convective T Tauri star BP Tau (0.7 Mq) 
exhibit a strong large-scale magnetic field (Donati ct al. 
2008a), whereas the more massive partly convective star 
V2129 Oph (1.4 Mq) possesses a weaker and more complex 
field (Donati ct al. 2007). This is reminiscent of the transi- 
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Figure 16. Magnetic flux as a function of Rossby number. On tlie left panel magnetic fluxes as measured from Stokes V spectra and 
ZDI by M08b and DOS arc mentioned as blue hexagons (M* > 0.4 Mq) and green squares (0.2 < M* < 0.4 Mq). Results from this paper 
(Aft Sj 0.2 Mq) are shown as red circles. On the right panel the ratio between the magnetic fluxes as recovered from Stokes V and Stokes / 
measurements (whenever available) are shown. Stokes / magnetic fluxes are taken from .lohns-Krull Sc Valcnti (2000), Rcincrs ,'(- Basri 
(2007), Rcincrs ct al. (2009), and Rcincrs & Basri (2009) (see Tab. 1). Measurements of different epochs (whenever available) for Stokes V 
are shown connected by a solid line. On each plot the Ro = 10~^, corresponding to the saturation level, is depicted as a vertical solid 
black line. Triangles denote upper limits. In the left panel, horizontal lines show the magnetic fluxes corresponding to saturation for stars 
with M* > 0.4 Mq (blue), and 0.2 < M* < 0.4 Mq (green). In the right panel they show the mean fraction of magnetic flux detected in 
Stokes V spectra. 



tion we observe at 0.5 Mq among main sequence M dwarfs. 
Recent observations, on the fully convective V2247 Oph 
(0.35 Mq) reveal a still weaker and more complex magnetic 
field (Douati ot al. 2010) which may correspond to the weak 
field late M dwarf we present here. 

We do not detect Stokes V signatures in individual spec- 
tra of the two latest stars of our sample VB 8 and VB 10. 
However the averaged LSD profile of VB 10 suggests the 
presence of a toroidal axisymmetric field component on this 
object. Further observations may confirm this first spectro- 
polarimetric detection on an ultra-cool dwarf. 

The first spectropolarimetric survey of M dwarfs has al- 
ready provided dynamo theorists with strong constraints on 
the evolution of surface magnetic fields of M dwarfs across 
the fully convective divide (DOS and M08b). The results pre- 
sented in this paper on the magnetic topologies of late M 
dwarfs reveal a new unexpected behaviour below 0.2 Mq. 
We interpret the fact that objects with similar stellar pa- 
rameters host radically different magnetic topologies as a 
possible evidence for a switch between two dynamo states 
(either cyclic or chaotic). Finally our observations suggest 
the presence of a large-scale magnetic field on the MS dwarf 
VB 10, featuring a significant toroidal axisymmetric com- 
ponent, whereas the detection of the magnetic field of VB 8 
(M7) is not possible from our spectropolarimetric data. 
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Table Al. Detailed journal of observations for GJ 51. Columns 1-7 list the UT date, the heliocentric Julian date, the UT time, the 
exposure time, the peak signal to noise ratio (per 2.6 kms~^ velocity bin) and the rms noise level (relative to the unpolarised continuum 
level and per 1.8 kms~^ velocity bin) in the average circular polarisation profile produced by Least-Squares Deconvolution (see Sec. 2.2). 
In column 7 we indicate the longitudinal field computed from Eq. 1. Column 8 lists the radial velocities (absolute accuracy 0.10 kms~^, 
internal accuracy 0.03 kms^^) associated to each exposure. The rotational cycle E from the cphemeris HJD = 2 453 950 + 1.02 E is 
mentioned in column 9. 



Date 


HJD 


UT 




t 


cxp 


S/N 


CLSD 


Bi 


RV 


Cycle 




(2,450,000+) 


(h:m:s) 




1 


;«) 




(10-4/c) 


(G) 


(kms-l) 




2006 






















Aug 05 


3953.13879 


15:15:58 


4 


X 


450.0 


121 


12.6 


-1617 ± 150 


-5.24 


3.077 


Aug 07 


3955.09204 


14:08:29 


4 


X 


550.0 


132 


9.3 


-1126 ± 100 


-5.37 


4.992 


Aug 08 


3956.08931 


14:04:27 


4 


X 


650.0 


165 


7.7 


-1105 ± 91 


-5.43 


5.970 


Aug 09 


3957.08253 


13:54:37 


4 


X 


500.0 


139 


9.0 


-782 ± 94 


-5.56 


6.944 


Aug 11 


3959.09703 


14:15:18 


4 


X 


500.0 


149 


8.2 


-673 ± 92 


-5.80 


8.919 


Aug 12 


3960.10164 


14:21:51 


4 


X 


500.0 


128 


10.1 


-702 ± 103 


-5.74 


9.904 


2007 






















Sep 28 


4371.97896 


11:22:59 


4 


X 


840.0 


171 


6.7 


-1394 ± 96 


-7.14 


413.705 


Sep 30 


4374.02264 


12:25:53 


4 


X 


840.0 


179 


6.6 


-1389 ± 94 


-7.22 


415.708 


Oct 01 


4374.94925 


10:40:12 


4 


X 


840.0 


163 


7.0 


-2005 ± 111 


-6.12 


416.617 


Oct 02 


4375.82543 


07:41:54 


4 


X 


840.0 


159 


7.0 


-1395 ± 101 


-5.57 


417.476 


Oct 02 


4375.94877 


10:39:31 


4 


X 


840.0 


174 


6.5 


-2011 ± 108 


-5.85 


417.597 


Oct 02 


4376.05346 


13:10:16 


4 


X 


840.0 


180 


6.5 


-1522 ± 97 


-7.37 


417.699 


Oct 03 


4376.82249 


07:37:40 


4 


X 


840.0 


163 


6.8 


-1429 ± 104 


-5.52 


418.453 


Oct 03 


4376.94931 


10:40:18 


4 


X 


900.0 


195 


5.7 


-2076 ± 108 


-5.55 


418.578 


Oct 03 


4377.05178 


13:07:51 


4 


X 


900.0 


198 


5.7 


-1693 ± 95 


-6.87 


418.678 


2008 






















Oct 14 


4753.82296 


07:38:32 


4 


X 


815.0 


170 


7.5 


-1372 ± 90 


-7.16 


788.062 


Oct 14 


4753.86728 


08:42:22 


4 


X 


815.0 


183 


6.9 


-1049 ± 79 


-7.29 


788.105 


Oct 14 


4753.96099 


10:57:19 


4 


X 


815.0 


157 


8.0 


-619 ± 82 


-6.63 


788.197 


Oct 14 


4754.05403 


13:11:17 


4 


X 


815.0 


159 


7.6 


-516 ± 81 


-6.24 


788.288 


Oct 17 


4756.85209 


08:20:35 


4 


X 


815.0 


119 


10.4 


-1525 ± 105 


-6.13 


791.031 


Oct 20 


4759.77630 


06:31:33 


4 


X 


815.0 


176 


7.0 


-1600 ± 84 


-5.63 


793.898 


Oct 20 


4759.85136 


08:19:38 


4 


X 


815.0 


159 


8.0 


-1699 ± 95 


-6.13 


793.972 


Oct 20 


4759.94401 


10:33:04 


4 


X 


815.0 


152 


8.3 


-1521 ± 93 


-6.84 


794.063 


Oct 20 


4760.03093 


12:38:13 


4 


X 


815.0 


167 


7.7 


-1072 ± 85 


-7.34 


794.148 
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Table A2. Detailed journal of observations for GJ 1156. See Table Al for more details. 
HJD = 2 453 950 + 0.491 E. 



Data are phased according to the ephemeris 



Date 


HJD 


UT 




Ccxp 


S/N 


CLSD 


Bi 


RV 


Cycle 




(2,450,000+) 


(h;m:s) 




(s) 




(10-*/c) 


(G) 


(kms-i) 




2007 




















Mar 02 


4162.05198 


13:08:00 


4 


X 600.0 


157 


8.5 


84 ± 67 


5.86 


431.878 


Mar 03 


4163.05912 


13:18:14 


4 


X 420.0 


120 


11.0 


177 ± 81 


5.86 


433.929 


Mar 04 


4164.04908 


13:03:43 


4 


X 700.0 


181 


6.8 


157 ± 59 


5.83 


435.945 


Mar 05 


4165.08932 


14:01:37 


4 


X 600.0 


165 


7.7 


-15 ± 65 


5.80 


438.064 


Mar 06 


4166.03032 


12:36:37 


4 


X 500.0 


143 


9.3 


91 ± 73 


5.89 


439.980 


Mar 07 


4167.00828 


12:04:50 


4 


X 600.0 


163 


7.9 


-2 ± 74 


6.50 


441.972 


2008 




















Jan 20 


4485.98764 


11:39:04 


4 


X 540.0 


158 


8.5 


-264 ± 63 


5.72 


1091.625 


Jan 20 


4486.14332 


15:23:14 


4 


X 540.0 


156 


8.4 


148 ± 61 


5.95 


1091.942 


Jan 21 


4486.98190 


11:30:41 


4 


X 540.0 


127 


10.4 


-224 ± 74 


5.61 


1093.649 


Jan 21 


4487.14392 


15:23:59 


4 


X 540.0 


151 


8.8 


26 ± 65 


5.95 


1093.979 


Jan 22 


4488.09159 


14:08:30 


4 


X 610.0 


150 


8.8 


80 ± 65 


5.83 


1095.910 


2009 




















Jan 08 


4840.01516 


12:20:05 


4 


X 815.0 


195 


6.5 


-187 ± 48 


5.78 


1812.658 


Jan 08 


4840.10476 


14:29:06 


4 


X 815.0 


187 


6.5 


55 ±48 


5.73 


1812.841 


Jan 08 


4840.17094 


16:04:23 


4 


X 815.0 


166 


8.8 


168 ± 67 


5.83 


1812.975 


Jan 09 


4841.03123 


12:43:05 


4 


X 815.0 


167 


7.7 


-66 ± 55 


5.10 


1814.728 


Jan 09 


4841.09873 


14:20:17 


4 


X 815.0 


188 


6.8 


44 ± 50 


5.77 


1814.865 


Jan 09 


4841.16680 


15:58:18 


4 


X 815.0 


184 


6.7 


88 ±49 


5.97 


1815.004 


Jan 10 


4842.02100 


12:28:15 


4 


X 815.0 


184 


6.8 


-105 ± 49 


5.76 


1816.743 


Jan 10 


4842.09687 


14:17:29 


4 


X 815.0 


183 


6.9 


81 ± 51 


5.73 


1816.898 


Jan 10 


4842.16403 


15:54:11 


4 


X 815.0 


187 


6.7 


135 ± 50 


5.94 


1817.035 



Table A3. Detailed journal of observations for GJ 1245 B. See Table Al for more details. Data are phased according to the ephemeris 
HJD = 2 453 950 + 0.71 £. 



Date 


HJD 


UT 




icxp 


S/N 


CLSD 


Bi 


RV 


Cycle 




(2,450,000+) 


(h:m:s) 




(s) 




(10-*/c) 


(G) 


(kms-i) 




2006 




















Aug 05 


3952.96705 


11:05:26 


4 


X 630.0 


166 


8.6 


-251 ± 33 


5.29 


4.179 


Aug 08 


3955.95886 


10:53:38 


4 


X 700.0 


191 


7.1 


-79 ± 26 


5.40 


8.393 


Aug 09 


3956.96600 


11:03:55 


4 


X 600.0 


175 


7.7 


186 ± 28 


5.56 


9.811 


Aug 10 


3957.95942 


10:54:27 


4 


X 700.0 


185 


7.5 


-179 ± 27 


5.28 


11.210 


Aug 11 


3958.96340 


11:00:11 


4 


X 700.0 


185 


7.5 


145 ± 27 


5.52 


12.625 


Aug 12 


3959.96518 


11:02:45 


4 


X 700.0 


158 


8.8 


-133 ± 29 


5.48 


14.035 


2007 




















Sep 30 


4373.89061 


09 


17:56 


4 


X 840.0 


197 


6.8 


-107 ± 27 


5.30 


597.029 


Oct 01 


4374.86550 


08 


41:51 


4 


X 840.0 


199 


6.6 


148 ± 27 


5.47 


598.402 


Oct 02 


4375.72600 


05 


21:04 


4 


X 840.0 


214 


6.4 


63 ± 26 


5.51 


599.614 


Oct 02 


4375.86847 


08 


46:14 


4 


X 840.0 


182 


7.4 


-136 ± 29 


5.26 


599.815 


Oct 03 


4376.72595 


05 


21:05 


4 


X 900.0 


226 


5.7 


-179 ± 25 


5.27 


601.022 


Oct 03 


4376.86630 


08 


43:13 


4 


X 900.0 


218 


6.0 


108 ± 26 


5.46 


601.220 


2008 




















Aug 20 


4698.75329 


05:57:45 


4 


X 715.0 


187 


7.3 


-44 ± 25 


5.38 


1054.582 


Aug 20 


4698.84202 


08:05:31 


4 


X 715.0 


191 


7.3 


-64 ± 26 


5.34 


1054.707 


Aug 20 


4698.99653 


11:48:01 


4 


X 715.0 


171 


8.2 


45 ± 28 


5.52 


1054.925 


Aug 21 


4699.75240 


05:56:29 


4 


X 715.0 


194 


7.2 


56 ± 25 


5.52 


1055.989 


Aug 21 


4699.87910 


08:58:56 


4 


X 715.0 


186 


7.2 


26 ± 26 


5.53 


1056.168 


Aug 21 


4699.99952 


11:52:21 


4 


X 715.0 


176 


8.0 


70 ± 27 


5.60 


1056.337 


Aug 22 


4700.74540 


05:46:26 


4 


X 715.0 


150 


9.3 


74 ± 32 


5.44 


1057.388 


Aug 22 


4700.87247 


08:49:25 


4 


X 715.0 


153 


9.3 


-138 ± 31 


5.29 


1057.567 


Aug 22 


4700.99137 


11:40:38 


4 


X 715.0 


151 


9.6 


-47 ± 32 


5.45 


1057.734 


Aug 22 


4701.02897 


12:34:46 


4 


X 715.0 


138 


10.1 


-41 ± 33 


5.48 


1057.787 
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Table A4. Detailed journal of observations for WX UMa. See Table Al for more details. Data arc phased according to the ephemeris 
HJD = 2 453 850 + 0.78£;. 



Date 


HJD 


UT 




t 


cxp 


S/N 


0"LSD 


Be 


RV 


Cycle 




(2,450,000+) 


(h:m;s) 






[s) 




(10-4/c) 


(G) 


(kms-i) 




2006 






















Jun 09 


3895.77128 


06 


31:12 


4 


X 


600.0 


86 


16.1 


-1128 ± 213 


70.39 


0.476 


Jun 09 


3895.83862 


08 


08:10 


4 


X 


600.0 


67 


19.8 


-914 ± 207 


70.60 


0.562 


Jun 10 


3896.77760 


06 


40:22 


4 


X 


600.0 


137 


11.9 


-1816 ± 162 


70.59 


1.766 


Jun 10 


3896.84666 


08 


19:49 


4 


X 


600.0 


120 


11.0 


-2150 ± 198 


70.78 


1.855 


Jun 11 


3897.76780 


06 


26:19 


4 


X 


600.0 


142 


9.6 


-2053 ± 174 


68.97 


3.036 


Jun 11 


3897.83460 


08 


02:31 


4 


X 


600.0 


130 


10.4 


-1731 ± 119 


70.21 


3.121 


Jun 12 


3898.74596 


05 


54:56 


4 


X 


600.0 


122 


11.2 


-1147 ± 101 


70.03 


4.290 


Jun 12 


3898.81286 


07 


31:16 


4 


X 


600.0 


120 


11.2 


-1110 ± 103 


70.43 


4.375 


2007 






















Mar 02 


4161.97902 


11:26:01 


4 


X 


700.0 


143 


8.8 


-1825 ± 116 


69.91 


341.768 


Mar 03 


4162.99705 


11:51:58 


4 


X 


600.0 


115 


11.8 


-1340 ± 124 


70.04 


343.073 


Mar 04 


4163.97675 


11:22:44 


4 


X 


700.0 


154 


8.4 


-1581 ± 105 


69.93 


344.329 


Mar 05 


4165.01491 


12:17:42 


4 


X 


700.0 


152 


8.4 


-2357 ± 124 


69.90 


345.660 


Mar 07 


4166.94370 


10:35:10 


4 


X 


560.0 


137 


10.0 


-1265 ± 107 


70.03 


348.133 


Mar 08 


4168.03986 


12:53:39 


4 


X 


600.0 


139 


9.3 


-2175 ± 125 


69.92 


349.538 


2008 






















Jan 19 


4484.92606 


10:10:43 


4 


X 


560.0 


63 


21.4 


-1416 ± 206 


70.56 


755.803 


Jan 20 


4486.11357 


14:40:40 


4 


X 


500.0 


116 


12.0 


-1725 ± 115 


70.01 


757.325 


Jan 21 


4486.92292 


10:06:06 


4 


X 


540.0 


123 


11.2 


-2135 ± 123 


69.97 


758.363 


Jan 21 


4487.11353 


14:40:34 


4 


X 


520.0 


129 


10.4 


-1968 ± 114 


70.05 


758.607 


2009 






















Jan 11 


4843.01564 


12:20:02 


4 


X 


715.0 


116 


11.9 


-1147 ± 108 


70.24 


1214.892 


Jan 11 


4843.11846 


14:48:06 


4 


X 


715.0 


113 


12.5 


-1151 ± 114 


70.02 


1215.024 


Jan 12 


4843.94502 


10:38:18 


4 


X 


715.0 


132 


10.0 


-1670 ± 100 


69.86 


1216.083 


Jan 12 


4844.04478 


13:01:57 


4 


X 


715.0 


113 


12.2 


-1707 ± 110 


70.04 


1216.211 


Jan 12 


4844.14627 


15:28:05 


4 


X 


715.0 


133 


10.0 


-1772 ± 98 


69.91 


1216.341 


Jan 13 


4844.94269 


10:34:54 


4 


X 


715.0 


133 


10.0 


-1638 ± 101 


69.87 


1217.362 


Jan 13 


4845.09115 


14:08:40 


4 


X 


715.0 


146 


9.0 


-1864 ± 101 


69.76 


1217.553 


Jan 13 


4845.16697 


15:57:51 


4 


X 


715.0 


153 


8.4 


-1735 ± 101 


69.62 


1217.650 


Jan 14 


4845.98095 


11:29:57 


4 


X 


715.0 


145 


9.2 


-1274 ± 90 


69.50 


1218.694 


Jan 14 


4846.08129 


13:54:25 


4 


X 


715.0 


163 


8.0 


-1366 ± 86 


69.57 


1218.822 


Jan 14 


4846.17078 


16:03:17 


4 


X 


715.0 


158 


8.4 


-1130 ± 89 


69.75 


1218.937 
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Table A5. Detailed journal of observations for DX Cnc. See Table Al for more details. 
HJD = 2 453 950 + 0.46 £. 



Data are phased according to the ephemeris 



Date 


HJD 


UT 




tcxp 


S/N 


CLSD 


Bi 


RV 


Cycle 




(2,450,000+) 


(h:m:s) 




(s) 




{10-*/c) 


(G) 


(kms-i) 




2007 




















Mar 02 


4161.94084 


10:30:16 


4 


X 700.0 


165 


8.8 


145 ± 49 


10.60 


460.741 


Mar 04 


4163.93885 


10:27:33 


4 


X 700.0 


170 


8.6 


158 ± 47 


10.52 


465.084 


Mar 05 


4164.94533 


10:36:57 


4 


X 700.0 


179 


8.3 


-13 ± 46 


10.65 


467.272 


Mar 06 


4165.93747 


10:25:43 


4 


X 560.0 


119 


12.3 


209 ± 65 


10.54 


469.429 


Mar 07 


4166.91257 


09:49:56 


4 


X 560.0 


139 


11.1 


161 ± 55 


10.46 


471.549 


2008 




















Jan 18 


4483.86547 


08:40:11 


4 


X 520.0 


116 


12.7 


81 ± 61 


10.85 


1160.577 


Jan 19 


4484.86323 


08:36:58 


4 


X 560.0 


141 


10.6 


34 ± 52 


10.88 


1162.746 


Jan 19 


4485.04885 


13:04:15 


4 


X 570.0 


90 


16.7 


86 ± 79 


11.00 


1163.150 


Jan 20 


4486.05051 


13:06:39 


4 


X 600.0 


161 


9.4 


11 ±46 


10.70 


1165.327 


Jan 21 


4486.86360 


08:37:30 


4 


X 560.0 


138 


10.9 


135 ± 54 


10.43 


1167.095 


Jan 21 


4487.04831 


13:03:28 


4 


X 700.0 


160 


9.8 


131 ± 58 


9.94 


1167.496 


Jan 22 


4487.98223 


11:28:19 


4 


X 610.0 


152 


9.5 


166 ± 54 


9.29 


1169.527 


2009 




















Feb 13 


4875.88471 


09:08:28 


4 


X 800.0 


122 


12.3 


86 ± 56 


10.76 


2012.793 


Feb 14 


4876.87737 


08:57:56 


4 


X 800.0 


137 


11.6 


149 ± 52 


10.80 


2014.951 


Feb 14 


4877.00055 


11:55:19 


4 


X 800.0 


106 


14.8 


-4± 65 


10.68 


2015.219 


Feb 15 


4877.75183 


05:57:12 


4 


X 800.0 


164 


9.2 


48 ±43 


10.66 


2016.852 


Feb 15 


4877.87789 


08:58:44 


4 


X 800.0 


180 


8.5 


13 ±40 


10.49 


2017.126 


Feb 15 


4878.01424 


12:15:05 


4 


X 800.0 


157 


9.4 


110 ± 44 


10.58 


2017.422 


Feb 16 


4878.76253 


06:12:40 


4 


X 800.0 


187 


7.8 


54 ± 37 


10.61 


2019.049 


Feb 16 


4878.88556 


09:09:50 


4 


X 800.0 


185 


7.9 


78 ± 36 


10.78 


2019.316 


Feb 16 


4879.00624 


12:03:37 


4 


X 800.0 


177 


8.4 


71 ± 38 


10.71 


2019.579 



Table A6. Detailed journal of observations for GJ 3622. See Table Al for more details. 
HJD = 2 453 950 + 1.5 B. 



Data are phased according to the ephemeris 



Date 


HJD 


UT 




t 


cxp 


S/N 


O-LSD 


Bi 


RV 


Cycle 




(2,450,000+) 


(h:m:s) 






(s) 




(10-4/e) 


(G) 


(kms-i) 




2008 






















Mar 20 


4545.85137 


08:17:36 


4 


X 


990.0 


159 


9.3 


-50 ± 17 


2.20 


397.234 


Mar 21 


4546.94418 


10:31:17 


4 


X 


990.0 


162 


9.1 


-2 ± 18 


1.77 


397.963 


Mar 22 


4548.00749 


12:02:30 


4 


X 


990.0 


154 


9.5 


-20 ± 18 


1.88 


398.672 


Mar 23 


4548.92606 


10:05:17 


4 


X 


990.0 


156 


9.4 


-72 ± 17 


2.33 


399.284 


Mar 25 


4550.82345 


07:37:37 


4 


X 


990.0 


152 


9.8 


-60 ± 18 


2.57 


400.549 


Mar 26 


4551.82107 


07:34:14 


4 


X 


990.0 


167 


8.9 


-56 ± 16 


2.53 


401.214 


Mar 27 


4552.91578 


09:50:41 


4 


X 


990.0 


150 


10.1 


7± 18 


2.46 


401.944 


Mar 28 


4553.86120 


08:32:08 


4 


X 


990.0 


128 


11.1 


-5 ± 20 


2.41 


402.574 


2009 






















Feb 13 


4875.93524 


10:18:39 


4 


X 


900.0 


101 


15.0 


-22 ± 26 


2.35 


617.290 


Feb 15 


4877.93049 


10:11:42 


4 


X 


900.0 


134 


11.0 


-73 ± 20 


2.44 


618.620 


Feb 16 


4878.93305 


10:15:21 


4 


X 


900.0 


151 


10.0 


-42 ± 18 


2.32 


619.289 


Feb 16 


4879.06535 


13:25:51 


4 


X 


900.0 


135 


11.3 


-15 ± 19 


2.37 


619.377 


Feb 17 


4879.91581 


09:50:29 


4 


X 


900.0 


162 


9.3 


13 ± 17 


2.32 


619.944 


Feb 17 


4880.05123 


13:05:29 


4 


X 


900.0 


151 


10.0 


-18 ± 18 


2.42 


620.034 
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Table A7. Detailed journal of observations for GJ 1154 A. See Table Al for more details. The rotational cycle is not indicated since no 
rotation period could be derived for this star. 



Date 


HJD 


UT 




Ccxp 


S/N 


O-LSD 


Be 


RV 




(2,450,000+) 


(h:m:s) 




(s) 




[lO-^Ic] 


(G) 


(kms-i) 


2007 


















Mar 02 


4162.01888 


12:19:34 


4 


X 600.0 


150 


7.9 


-825 ± 63 


-12.76 


Mar 03 


4163.03213 


12:38:35 


4 


X 500.0 


118 


10.5 


-672 ± 69 


-12.79 


Mar 04 


4164.01521 


12:14:11 


4 


X 600.0 


167 


7.0 


-718 ± 55 


-12.77 


Mar 05 


4165.05486 


13:11:23 


4 


X 600.0 


147 


7.2 


-612 ± 55 


-13.07 


Mar 06 


4166.00262 


11:55:57 


4 


X 500.0 


152 


8.1 


-660 ± 59 


-12.82 


Mar 07 


4166.97810 


11:20:36 


4 


X 500.0 


140 


8.6 


-800 ± 62 


-12.78 


2008 


















Jan 19 


4485.01555 


12:18:53 


4 


X 540.0 


86 


15.9 


-576 ± 85 


-12.81 


Jan 20 


4486.01710 


12:20:59 


4 


X 500.0 


147 


8.2 


-776 ± 57 


-12.89 


Jan 21 


4487.01179 


12:13:12 


4 


X 540.0 


130 


9.3 


-721 ± 61 


-12.76 


Jan 22 


4488.01494 


12:17:37 


4 


X 610.0 


154 


8.3 


-727 ± 61 


-13.22 



Table A8. Detailed journal of observations for GJ 1224. See Table Al for more details. The rotational cycle is not indicated since no 
rotation period could be derived for this star. 



Date 


HJD 


UT 




Ccxp 


S/N 


O-LSD 


Bi 


RV 




(2,450,000+) 


(h:m:s) 




(s) 




(10-4/,) 


(G) 


(kms-i) 


2008 


















Jun 25 


4643.02829 


12:33:21 


4 


X 800.0 


176 


6.2 


-535 ± 36 


-32.75 


Jun 28 


4645.88511 


09:07:15 


4 


X 800.0 


51 


19.7 


-546 ± 95 


-32.69 


Jun 28 


4645.92645 


10:06:46 


4 


X 800.0 


86 


14.1 


-631 ± 67 


-32.62 


Jun 29 


4646.94450 


10:32:47 


4 


X 800.0 


147 


7.9 


-550 ± 42 


-32.67 


Jul 02 


4649.94447 


10:32:50 


4 


X 800.0 


155 


7.3 


-632 ± 41 


-32.72 


Jul 03 


4650.91640 


09:52:27 


4 


X 800.0 


171 


6.6 


-550 ± 37 


-32.69 


Jul 04 


4651.85383 


08:22:23 


4 


X 800.0 


171 


6.3 


-556 ± 37 


-32.64 


Jul 05 


4652.86038 


08:31:51 


4 


X 800.0 


176 


6.1 


-617 ± 37 


-32.69 


Jul 06 


4653.83070 


07:49:09 


4 


X 800.0 


185 


5.8 


-545 ± 34 


-32.64 


Jul 07 


4654.89004 


09:14:39 


4 


X 800.0 


166 


6.6 


-573 ± 38 


-32.65 


Jul 08 


4655.84785 


08:13:57 


4 


X 800.0 


172 


6.4 


-523 ± 37 


-32.63 


Jul 09 


4656.96336 


11:00:20 


4 


X 800.0 


172 


6.4 


-500 ± 40 


-32.72 



Table A9. Detailed journal of observations for CN Leo. See Table Al for more details. The rotational cycle is not indicated since no 
rotation period could be derived for this star. 



Date 



HJD 

(2,450,000^ 



UT 
(h:m:s) 



(s) 



S/N (TLSD 

(10-'4/c 



Bi RV 

(G) (kms~i) 



2008 
















Mar 21 


4546.99345 


11:43:06 


4 X 800.0 


179 


5.6 


-651 ± 45 


19.66 


Mar 22 


4547.91616 


09:51:50 


4 X 800.0 


172 


6.7 


-635 ± 48 


19.68 


Mar 22 


4547.95840 


10:52:40 


4 X 800.0 


186 


6.1 


-708 ± 46 


19.63 


Mar 23 


4548.87485 


08:52:23 


4 X 800.0 


245 


4.5 


-772 ± 45 


19.55 
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Table AlO. Detailed journal of observations for VB 8. See Table Al for more details. The rotational cycle is not indicated since no 
rotation period could be derived for this star. 



Date 


HJD 


UT 




Ccxp 


S/N 


O-LSD 


Bi 


RV 




(2,450,000+) 


(h:m:s) 




(s) 




(lO-^/c) 


(G) 


(kms-i) 


2009 


















May 03 


4954.89250 


09:18:10 


4 


X 1154.0 


104 


16.5 


24 ± 29 


15.41 


May 03 


4954.99885 


11:51:19 


4 


X 1154.0 


107 


15.8 


-15 ± 28 


15.24 


May 03 


4955.10357 


14:22:06 


4 


X 1154.0 


107 


15.7 


-52 ± 28 


15.39 


May 04 


4955.90408 


09:34:47 


4 


X 1154.0 


83 


20.0 


113 ± 36 


15.63 


May 04 


4956.00991 


12:07:11 


4 


X 1154.0 


94 


17.9 


49 ± 32 


15.41 


May 04 


4956.11013 


14:31:29 


4 


X 1154.0 


98 


17.0 


-25 ± 30 


15.31 


May 05 


4956.90261 


09:32:37 


4 


X 1154.0 


99 


16.9 


77 ± 30 


15.31 


May 05 


4957.00469 


11:59:36 


4 


X 1154.0 


102 


16.2 


87 ± 29 


15.32 


May 05 


4957.11128 


14:33:05 


4 


X 1154.0 


103 


16.3 


5 ± 29 


15.51 



Table All. Detailed journal of observations for VB 10. See Table Al for more details. The rotational cycle is given for the most probable 
period according to the ephemeris HJD = 2 455 000 -I- 0.69 E. 



Date 


HJD 


UT 




Ccxp 


S/N 


0"LSD 


Bi 


RV 


Cycle 




(2,450,000+) 


(h:m:s) 




(s) 




(10-*7c) 


(G) 


(kms-i) 




2009 




















Jul 01 


5013.96290 


10:58:52 


4 


X 1154.0 


68 


25.7 


94 ± 65 


36.16 


20.236 


Jul 01 


5014.01972 


12:20:41 


4 


X 1154.0 


77 


23.0 


92 + 57 


36.25 


20.318 


Jul 01 


5014.07652 


13:42:29 


4 


X 1154.0 


76 


23.0 


145 ± 57 


36.03 


20.401 


Jul 02 


5014.95765 


10:51:17 


4 


X 1154.0 


71 


24.2 


16 + 59 


36.25 


21.678 


Jul 02 


5015.01489 


12:13:42 


4 


X 1154.0 


74 


23.4 


68 ± 57 


36.37 


21.761 


Jul 02 


5015.10842 


14:28:23 


4 


X 1154.0 


80 


22.5 


93 ± 58 


36.13 


21.896 


Jul 03 


5015.94279 


10:29:52 


4 


X 1154.0 


74 


23.1 


-69 ± 59 


36.27 


23.105 


Jul 03 


5016.04819 


13:01:38 


4 


X 1154.0 


75 


23.6 


-4 + 56 


36.51 


23.258 


Jul 03 


5016.10641 


14:25:28 


4 


X 1154.0 


74 


23.5 


92 ± 61 


36.13 


23.343 



